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The flavor composition of high-energy astrophysical neutrinos can reveal the physics governing 
their production, propagation, and interaction. The IceCube Collaboration has published the first 
experimental determination of the ratio of the flux in each flavor to the total. We present, as a 
theoretical counterpart, new results for the allowed ranges of flavor ratios at Earth for arbitrary 
flavor ratios in the sources. Our results will allow IceCube to more quickly identify when their data 
imply standard physics, a general class of new physics with arbitrary (incoherent) combinations 
of mass eigenstates, or new physics that goes beyond that, e.g., with terms that dominate the 
Hamiltonian at high energy. 


Introduction. — The discovery of astrophysical neu¬ 
trinos with energies up to a few PeV by the IceCube 
Collaboration is tremendously important for multi¬ 
messenger astronomy as well as for new tests of neutrino 
properties. While the origin of these neutrinos is still un¬ 
clear, there are important clues in the energy spectrum 
and sky distribution, and a component from cosmic dis¬ 
tances (^ Gpc) is required [SH2l- These are the most 
extreme energies and distances for detected neutrinos. 

The flavor composition is also expected to be impor¬ 
tant, because the ratio of flux in each flavor to the to¬ 
tal cancels the unknown normalization. The ratios de¬ 
pend on the physical conditions at the source, the effects 
of standard flavor mixing, and on potential new physics 

I51E5H361. 

The first IceCube results on flavor composition have 
been published recently m, and were followed by results 
obtained with a combined-likelihood analysis of several 
data sets with more statistics [37]. Accordingly, there 
has been intense interest in deducing flavor ratios from 
IceCube data |9l|3ll|3l|3i[39]. 

In this Letter, we use ternary plots or “flavor triangles” 
to show the flavor composition at Earth. We system¬ 
atically explore which regions of this plot can be pop¬ 
ulated from theoretical perspectives —without or with 
new physics— including the uncertainties in source fla¬ 
vor composition and neutrino mixing parameters. We 
also note prospects for the proposed volume upgrade, 
IceCube-Gen 2 [40] . 

We make no distinction between u and P, because, 
except for yet-unobserved high-energy events, IceCube 
cannot distinguish between them. (In addition, their 
cross sections agree to better than 5% in this energy 
range [1TJ|42].) 

All plots shown in the main text are for the normal 
neutrino mass hierarchy (NH), in which ui is the lightest 
mass eigenstate. Corresponding plots for the inverted 
hierarchy (IH), in which z /3 is lightest, are given in the 
Supplemental Material; the differences are modest. 


Flavor identification in IceCube. — IceCube can 
discriminate between muon tracks (from z/^, mostly) and 
cascades (from charged-current interactions of z/g and z/^-, 
mainly, and from neutral-current interactions of all fla¬ 
vors). If higher-energy events are observed, it will be 
possible to isolate i>e cascades via the Glashow reso¬ 
nance [43II45] . and z/^- and via double-bang and lol¬ 
lipop topologies [46H48] . In their absence, there is an 
experimental degeneracy between the electron and tau 
neutrino flavor content at Earth [34l [35]. In contrast, 
theoretically predicted flavor ratios, even in models with 
new physics, have a /i-r symmetry due to that mixing 
angle being near-maximal. 

Flavor composition at the source. — The flavor 
composition at the source could be quite different de¬ 
pending on the physical conditions. For the pion decay 
chain, which is often used as standard (“pion beam”), one 
expects a composition (/e,s : f^l,s ■ /t,s)=(| : § : 0)s, 
with fa,s the ratio of to the total flux, where 

/e,s + f^,s + /r,s = 1- Synchrotron cooling of secondary 
muons in strong magnetic fields leads to a transition to 
(0 : 1 : 0)s (“muon damped”) at higher energies, which 
depends on the field strength; see, e.g., Refs. 0131 Ei- 
im. If these muons pile up at lower energies [5l], or 
if there are contributions from charmed meson decays 
[29[|52j|53], then ^ ' 0)s is expected. Neutron de¬ 
cays [ 5 ] lead to (1 : 0 : 0)s. Small deviations, < 5% in 
the Vejv^ ratio, are expected from effects such as the he- 
licity dependence of muon decays 0IS3]. If several of the 
above processes in the source compete, arbitrary flavor 
compositions (/e,s • 1 — /e,s • 0) can be obtained [51]. If, 
in addition, are produced, such as by oscillations in a 
matter envelope even (/e,s : /^,s : 1 -/e,s -/ai.s) 

(with 0 < < 1 — /e,s) could be possible. Dark mat¬ 

ter annihilation or decay could yield any mixture, but 
(^ : ^ : |)s is the most natural. 

Flavor composition at Earth. — Here we focus 
on a diffuse flux, which is composed of small contribu¬ 
tions from many sources over a wide range of distances. 
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FIG. 1. Flavor content of the three active mass eigenstates. 
The regions are given by the best-fit values of the mixing pa¬ 
rameters (light yellow), and their Icr (darker) and Scr (dark¬ 
est) uncertainty regions [62], assuming a normal mass hierar¬ 
chy (NH). The tilt of the tick marks indicates the orientation 
with which to read the flavor content. 

and detected with energy resolution > 10 % (and binned 
more coarsely). In this case, the neutrinos are, at least 
effectively, an incoherent mixture of mass eigenstates. 
Even for the solar AttIq « 8 • 10“^ eV^ and PeV ener¬ 
gies, the vacuum oscillation length is only ^ 10 “^^ Gpc, 
much smaller than the complete baseline. (Depending 
on the physics in the production region, there can be 
also wave packet decoherence in the source [58f[6Q] .) As 
a consequence, the flavor composition at Earth [58] is 
// 3 ,e = Ei a /a.s> with U the PMNS matrix 

m, implying = 1. Eor a pion beam, the flavor 

composition evolves roughly into flavor equipartition at 
the detector, (| : | : |) 0 . 

New physics in neutrino propagation might modify the 
flavor composition. We categorize classes of new-physics 
models below. 

Flavor content of the mass eigenstates.— 

Eigure shows the flavor content \Uai\‘^ of the mass 
eigenstates, which is the fundamental input that deter¬ 
mines flavor ratios at Earth without or with new physics. 
It also illustrates the underlying three-flavor unitarity of 
our analysis, z.e., \Uai\‘^ + which 

allows the flavor content to be displayed in a ternary 
plot [63]. This is appropriate because the mixing angles 
to sterile neutrinos must be quite small [Ml ES] • 

The long axis of each region is set by the uncertainty in 
6^23 and (5cp, while the short axis is set by the uncertainty 
in 6 ^ 12 . The effect of the uncertainty in Ois is tiny. Even if 
^23 were to be precisely determined soon, it is less likely 
that ^cp will be, and the uncertainty in the latter will 


FIG. 2. Allowed flavor ratios at Earth with no new physics. 
The flavor ratios at the source are arbitrary (gray) or contain 
no tan flavor (red). The IceGube results are from Ref. [37] . 

still span a large range in \Urif and \Ur 2 f‘ 

Standard flavor mixing. — Eigure shows the al¬ 
lowed region for the flavor composition at Earth assum¬ 
ing arbitrary flavor composition at the source and stan¬ 
dard neutrino mixing (including parameter uncertain¬ 
ties). The region is quite small: even at 3cr it covers 
only about 10% of the available space. There is little 
difference between /r,s = 0 and fr,s % 0. 

There is a theoretical symmetry along the line (/e,© • 
(1 — / e ,©)/2 : (1 — / e ,©)/ 2 ) coming from nearly-maximal 
mixing. On the other hand, the experimental degener- 
acy pulls towards (/e,e : : 1 - - fe,®), with 

/e,© < 1 — /^,©, on account of the difficulty of distin¬ 
guishing between electromagnetic and hadronic cascades. 
Thus, theory and experiment are complementary, which 
enhances the discriminating power of flavor ratios. 

The region shown includes the possibility of energy- 
dependent flavor composition at the source; see the Sup¬ 
plemental Material for an example. It also includes the 
possibility that the diffuse flux has contributions from 
sources with different flavor compositions, because of the 
linear mapping between those at the source and those at 
Earth. 

Whereas the first IceCube flavor ratio analysis [35] 
used only three years of contained-vertex events, the up¬ 
dated analysis whose exclusion curves are shown in 
Fig.[2l combines several different data sets collected over 
four years, including through-going muons. The exclu¬ 
sion curves of both analyses are compatible. 

Eigure shows that if the flavor composition at the 
source could be restricted from astrophysical arguments, 
the allowed regions at Earth could become tiny (and will 
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FIG. 3. Allowed flavor ratios at Earth for different choices of 
source ratios, assuming standard mixing. Projected Icr, 2cr, 
and 3(7 exclusion curves from IceCube-Gen2 are included for 
comparison (gray, dotted); see main text. 


shrink when the mixing parameters are better known). A 
source composition of (1 : 0 : 0 )g is already disfavored at 
> 2 ( 7 . While the current IceCube fit is compatible with 
the standard | at la, the best-fit point cannot 

be reached within the Standard Model. 

An upgrade of IceCube would have excellent discrim¬ 
ination power, as indicated by the projected sensitivity 
curves we estimate for IceCube-Gen 2 and show in Fig.|^ 
We reduced the IceCube uncertainties by a factor 5, cor¬ 
responding to an exposure increased by a factor ^ 25 
6 times larger effective area m and twelve years 
instead of three). The true sensitivity might be worse 
(due to sparser instrumentation) or better (due to new 
techniques or to the discovery of flavor-identifying sig¬ 
nals [44l [46l HSl [SU [66H74) ). To be conservative, 
we assumed the best fit will correspond to the most- 
frequently considered composition, (| : ^ : |) 0 , for 
which it will be most difficult to test for new physics. 

Flavor ratios with new physics. — New physics 
can modify the flavor composition at production, during 
propagation, or in interaction. In the first two cases, it 
will affect the flavor composition that reaches the detec¬ 
tor; this is our focus. In the last case —which includes, 
e.g.^ non-standard interactions UB and renormalization 
group running of the mixing parameters m — we as¬ 
sume that new physics, possibly energy-dependent, can 
be separated by probing the interaction length in Earth 
via the angular dependence of the neutrino flux [TThSQ] . 

In extreme scenarios, there could be only one mass 
eigenstate present at detection, and the flavor composi¬ 
tion would correspond to that of one eigenstate. This 


EIG. 4. Allowed flavor ratios at Earth in a general class of 
new-physics models. These produce linear combinations of 
the flavor content of 1 ^ 3 , 1 ^ 2 , and z^i, shown as yellow (dashed) 
curves, from left to right. The standard mixing 3a region 
from Eig.[^is shown as a magenta (dotted) curve. 


could happen if all but one mass eigenstate completely 
decays or if matter-affected mixing at the source singles 
out a specific one for emission. 

Figure shows the allowed region if we restrict our¬ 
selves to a general class of new-physics models —those in 
which arbitrary combinations of incoherent mass eigen¬ 
states are allowed (we give examples below of mod¬ 
els that can access the area outside this region). The 
(a-flavor content of an allowed point is computed as 
\Uaif + ^2 1 ^ 4 ^ 21 ^ + ^3 ^ whcrc the ki are varied 

under the constraint + /c 2 + ^3 = 1 and the values of the 
mixing parameters are fixed. To generate the complete 
region, we repeat the procedure by varying the mixing 
parameters within their uncertainties. 

For a particular new-physics model, the functional 
forms and values of the ki are determined by its param¬ 
eters. The most dramatic examples include all variants 
of neutrino decay among mass eigenstates, both partial 
and complete [25l ISTHSd] . and secret neutrino interac¬ 
tions [HSHST]; the ki in these cases depend on neutrino 
lifetimes and new coupling constants, respectively. Other 
examples are pseudo-Dirac neutrinos j92f[9^ and deco¬ 
herence on the Planck-scale structure of spacetime [95]- 

[TOT] . 

Even with this general class of new-physics models, 
only about 25% of the flavor triangle can be accessed. 
The current IceCube best fit cannot be reached even by 
invoking this class of physics models. IceCube-Gen2 will 
be needed to strongly constrain such new-physics models. 

Interestingly, there is more than one way in which 






















4 


the standard ^ composition can be generated, 

such as through the standard mixing of (| : | : 0)g, or 
through a fortuitous incoherent mix of mass eigenstates 
due to decay. 

Already, complete decay in the most often used neu¬ 
trino decay scenario (only ui stable) for the NH can be 
ruled out at > 2cr (see Ref. [84] for a weaker exclusion at 
Icr based on their own analysis of tracks and cascades), 
and bounds on the neutrino lifetimes can be set m- 

To access the white region in Fig. [^ a broader class 
of new-physics models is required. Possible examples are 
models with violation of CPT and/or Lorentz invariance 
(which alter the dispersion relations) [25] llQQf llQ3HlQ6j , 
or the equivalence principle [10714109] . and coupling to a 
torsion field mni. 

All these have in common that they either invalidate 
the concept of decoherence in the astrophysical neutrino 
flavor composition or they change the values of the mix¬ 
ing parameters. Ref. m adopted a generic effective 
theory approach in which the new-physics terms domi¬ 
nate the propagation Hamiltonian at high energies, and 
showed that such models are indeed able to populate al¬ 
most the full triangle. 

Another possibility is the existence of extra dimen¬ 
sions, which could lead to matter-like resonant mixing 
between active and sterile flavors ma. Boosted dark 
matter [Blmlini] could generate neutrino-like events, 
even mimicking pure-flavor signatures. 

Conclusions. — We have demonstrated that the al¬ 
lowed region of neutrino flavor composition at Earth un¬ 
der standard mixing is quite small, in spite of the uncer¬ 
tainties in the mixing parameters and flavor composition 
at the sources. The allowed region remains small even 
in the presence of a general class of new-physics mod¬ 
els whose effect is to change the incoherent mix of mass 
eigenstates during propagation (e.^., neutrino decay and 
secret interactions). These results hardly depend on the 
mass hierarchy, and they hold for energy-dependent fla¬ 
vor compositions at the source or energy-dependent new 
physics, even when simultaneously present Ea; see the 
Supplemental Material. 

In order to access the larger space of possible fla¬ 
vor combinations, a broader class of new physics dur¬ 
ing propagation —flavor-violating or capable of modify¬ 
ing the values of the mixing parameters— or at detec¬ 
tion is required. Interestingly, the current IceCube best- 
fit composition lies in this region, though the standard 
(^ : I : 1)0 case is not excluded. 

The power of IceCube to determine the composition 
is enhanced by the complementarity between its experi¬ 
mental z^e-^r degeneracy and the theoretical sym¬ 

metry coming from nearly-maximal mixing. The current 
bounds are not only compatible with most source compo¬ 
sitions, but also with many potential new physics effects. 
However, the most favored neutrino decay scenario (only 
ui stable) can be already ruled out at > 2cr. 

The smaller the allowed region with only standard mix¬ 
ing shown in Fig. I^and Fig-i the more sensitive IceCube 


is to new physics. Likewise, the smaller the new-physics 
region shown in Fig. the more sensitive IceCube is to 
the broader class of new physics. The recent successes in 
measuring neutrino mixing parameters have been essen¬ 
tial to making these regions small. Our results provide 
new perspectives that will sharpen and accelerate tests 
of flavor ratios. 

Ideally, flavor ratios would be determined using a single 
class of point sources at known distances. No high-energy 
astrophysical sources have been resolved yet, however. 
We have shown that, even using a diffuse flux, flavor 
ratios can reveal information about source conditions and 
neutrino properties. 

Data from a volume upgrade of IceCube in combi¬ 
nation with improved measurements of the mixing pa¬ 
rameters, including (5cp, have the potential to nail down 
the flavor composition at the source or to identify new 
physics in propagation. However, it is not possible to 
extract the value of 6 cp from astrophysical data alone if 
the flavor composition at the source is not known; see 
the Supplemental Material. 

To fully exploit the power of neutrino flavors, advances 
in four directions are needed: 

1. A volume upgrade of IceCube (IceCube-Gen2 [40]) 
or a corresponding experiment in seawater (e.^., 
KMSNeT [TT6] L 

2. Reduction of the uncertainties in the values of the 
mixing parameters (especially O 23 and (5cp). 

3. Improvements in experimental techniques to recon¬ 
struct neutrino flavor and energy. 

4. More systematic model building to better under¬ 
stand, or constrain, the region of flavor ratios at 
Earth that could be accessed by new physics. 

Given the wealth of information about neutrino pro¬ 
duction, propagation, and interaction that the flavor 
composition provides, its precise determination should 
become a high-priority goal of ongoing and near-future 
experimental analyses. 

Acknowledgements. We thank Markus Ackermann, 
Shunsaku Horiuchi, Kohta Murase, Kenny Ng, Nathan 
Whitehorn, and Guanying Zhu for useful discussions and 
comments; we thank Marek Kowalski, Lars Mohr maun 
for that and for help reproducing the curves in Ref. m- 
JFB is supported by NSF Grant PHY-1404311. MB 
thanks the Institute for Nuclear Theory at the University 
of Washington for its hospitality during 14 - 26 June, 
2015, and the Department of Energy for partial sup¬ 
port during the completion of this work. WW thanks 
CCAPP for the hospitality during 29 April - 3 May, 2015, 
where substantial progress was made on this project, 
and acknowledges support from the “Helmholtz Alliance 
for Astroparticle Physics HAP” funded by the Initia¬ 
tive and Networking Fund of the Helmholtz Association. 
This project has received funding from the European 









Research Council (ERC) under the European Union’s no. 646623). 
Horizon 2020 research and innovation programme (grant 


5 


[1] M.G. Aartsen et al. (IceCube Collaboration), “First 
observation of PeV-energy neutrinos with IceCube,” 
Phys. Rev. Lett. Ill, 021103 (2013), arXiv: 1304.5356 
[astro-ph.HE] 

[2] M.C. Aartsen et al (IceCube), “Evidence for High- 
Energy Extraterrestrial Neutrinos at the IceCube De¬ 
tector,” Science 342 , 1242856 (2013), arXiv:1311.5238 
[astro-ph.HE] 

[3] M.G. Aartsen et al. (IceCube), “Search for a diffuse flux 
of astrophysical muon neutrinos with the IceCube 59- 
string configuration,” Phys. Rev. D89, 062007 (2014), 
arXiv:1311.7048 [astro-ph.HE] 

[4] M.G. Aartsen et al. (IceCube), “Observation of High- 
Energy Astrophysical Neutrinos in Three Years of Ice- 
Cube Data,” Phys. Rev. Lett. 113 , 101101 (2014), 
arXiv: 1405.5303 [astro-ph.HE] 

[5] Paolo Lipari, Maurizio Lusignoli, and Davide Melon!, 
“Flavor Composition and Energy Spectrum of Astro- 
physical Neutrinos,” Phys. Rev. D75, 123005 (2007), 
arXiv:0704.0718 [astro-ph] 

[6] Kohta Murase and Kunihito loka, “TeV-PeV Neutrinos 
from Low-Power Gamma-Ray Burst Jets inside Stars,” 
Phys. Rev. Lett. Ill, 121102 (2013), arXiv: 1306.2274 
[astro-ph.HE] 

[7] Ranjan Laha, John F. Beacom, Basudeb Dasgupta, 
Shunsaku Horiuchi, and Kohta Murase, “Demysti¬ 
fying the PeV Cascades in IceCube: Less (Energy) 
is More (Events),” Phys. Rev. D88, 043009 (2013), 
arXiv: 1306.2309 [astro-ph.HE] 

[8] Kohta Murase, Markus Ahlers, and Brian C. Lacki, 
“Testing the Hadronuclear Origin of PeV Neutrinos Ob¬ 
served with IceCube,” Phys. Rev. D88, 121301 (2013), 
arXiv: 1306.3417 [astro-ph.HE] 

[9] Walter Winter, “Photohadronic Origin of the TeV- 
PeV Neutrinos Observed in IceCube,” Phys. Rev. D88, 
083007 (2013), arXiv: 1307.2793 [astro-ph.HE] 

[10] Arman Esmaili and Pasquale Dario Serpico, “Are Ice- 
Cube neutrinos unveiling PeV-scale decaying dark mat¬ 
ter?” JCAP 1311 , 054 (2013), arXiv:1308.1105 [hep- 

[11] Soebur Razzaque, “The Galactic Center Origin of a 
Subset of IceCube Neutrino Events,” Phys. Rev. D88, 
081302 (2013), arXiv:1309.2756 [astro-ph.HE] 

[12] Markus Ahlers and Kohta Murase, “Probing the Galac¬ 
tic Origin of the IceCube Excess with Gamma-Rays,” 
Phys. Rev. D90, 023010 (2014), arXiv:1309.4077 [astro- 
ph.HE] 

[13] Ruo-Yu Liu, Xiang-Yu Wang, Susumu Inoue, Roland 
Crocker, and Felix Aharonian, “Diffuse PeV neutri¬ 
nos from EeV cosmic ray sources: Semirelativistic hy¬ 
pernova remnants in star-forming galaxies,” Phys. Rev. 
D89, 083004 (2014), arXiv:1310.1263 [astro-ph.HE] 

[14] Cecilia Lunardini, Soebur Razzaque, Kristopher T. 
Theodoseau, and Lili Yang, “Neutrino Events at Ice- 
Cube and the Fermi Bubbles,” Phys. Rev. D90, 023016 
(2014), arXiv:1311.7188 [astro-ph.HE] 

[15] Kazumi Kashiyama and Peter Meszaros, “Galaxy Merg¬ 


ers as a Source of Cosmic Rays, Neutrinos, and Gamma 
Rays,” Astrophys. J. 790, L14 (2014), arXiv:1405.3262 
[astro-ph.HE] 

[16] Xiao-Chuan Chang and Xiang-Yu Wang, “The diffuse 
gamma-ray flux associated with sub-PeV/PeV neutri¬ 
nos from starburst galaxies,” Astrophys. J. 793, 131 
(2014), arXiv: 1406.1099 [astro-ph.HE] 

[17] M.G. Aartsen et al. (IceCube), “Searches for Ex¬ 
tended and Point-like Neutrino Sources with Four Years 
of IceCube Data,” Astrophys. J. 796, 109 (2014), 
arXiv: 1406.6757 [astro-ph.HE] 

[18] Atri Bhattacharya, Rikard Enberg, Mary Hall Reno, 
and Ina Sarcevic, “Charm decay in slow-jet supernovae 
as the origin of the IceCube ultra-high energy neu¬ 
trino events,” JCAP 1506, 034 (2015), arXiv:1407.2985 
[astro-ph.HE] 

[19] Atri Bhattacharya, Raj Gandhi, and Aritra Gupta, 
“The Direct Detection of Boosted Dark Matter at High 
Energies and PeV events at IceCube,” JCAP 1503, 027 
(2015), arXiv: 1407.3280 [hep-ph] 

[20] M.G. Aartsen et al. (IceCube), “Searches for small- 
scale anisotropies from neutrino point sources with three 
years of IceCube data,” Astropart. Phys. 66, 39-52 
(2015), arXiv: 1408.0634 [astro-ph.HE] 

[21] Kohta Murase, “On the Origin of High-Energy Cos¬ 
mic Neutrinos,” Proeeedings of the 26th International 
Conferenee on Neutrino Physies and Astrophysies (Neu¬ 
trino 2014 ), AIP Conf. Proc. 1666, 040006 (2015), 
arXiv: 1410.3680 [hep-ph] 

[22] Luis A. Anchordoqui, “Neutron /3-decay as the origin 
of IceCubes PeV (anti)neutrinos,” Phys. Rev. D91, 
027301 (2015), arXiv: 1411.6457 [astro-ph.HE] 

[23] M.G. Aartsen et al. (IceCube), “Search for Prompt Neu¬ 
trino Emission from Gamma-Ray Bursts with IceCube,” 
Astrophys. J. 805, L5 (2015), arXiv:1412.6510 [astro- 
ph.HE] 

[24] Markus Ahlers, Yang Bai, Vernon Barger, and 
Ran Lu, “Galactic TeV-PeV Neutrinos,” (2015), 
arXiv:1505.03156 [hep-ph] 

[25] Gabriela Barenboim and Chris Quigg, “Neutrino ob¬ 
servatories can characterize cosmic sources and neu¬ 
trino properties,” Phys. Rev. D67, 073024 (2003), hep- 
ph/0301220 

[26] Zhi-zhong Xing and Shun Zhou, “Towards determina¬ 
tion of the initial flavor composition of ultrahigh-energy 
neutrino fluxes with neutrino telescopes,” Phys. Rev. 
D74, 013010 (2006), astro-ph/0603781 

[27] Sandip Pakvasa, Werner Rodejohann, and Thomas J. 
Weiler, “Flavor Ratios of Astrophysical Neutrinos: Im¬ 
plications for Precision Measurements,” JHEP 02, 005 
(2008), arXiv:0711.4517 [hep-ph] 

[28] Arman Esmaili and Yasaman Farzan, “An Analysis of 
Cosmic Neutrinos: Flavor Composition at Source and 
Neutrino Mixing Parameters,” Nucl. Phys. B821, 197- 
214 (2009), arXiv:0905.0259 [hep-ph] 

[29] Sandhya Choubey and Werner Rodejohann, “Flavor 
Composition of UHE Neutrinos at Source and at Neu- 







6 


trino Telescopes,” Phys. Rev. D80, 113006 (2009), 
arXiv:0909.1219 [hep-ph] 

[30] Kwang-Chang Lai, Guey-Lin Lin, and T.C. Liu, “De¬ 
termination of the Neutrino Flavor Ratio at the As- 
trophysical Source,” Phys. Rev. D80, 103005 (2009), 
arXiv:0905.4003 [hep-ph] 

[31] Olga Mena, Sergio Palomares-Ruiz, and Aaron C. Vin¬ 
cent, “Flavor Composition of the High-Energy Neu¬ 
trino Events in IceCube,” Phys. Rev. Lett. 113, 091103 
(2014), arXiv:1404.0017 [astro-ph.HE] 

[32] Xun-Jie Xu, Hong-Jian He, and Werner Rodejohann, 
“Constraining Astrophysical Neutrino Elavor Composi¬ 
tion from Leptonic Unitarity,” JCAP 1412, 039 (2014), 
arXiv: 1407.3736 [hep-ph] 

[33] Lingjun Eu, Chiu Man Ho, and Thomas J. Weiler, 
“Aspects of the Elavor Triangle for Cosmic Neu¬ 
trino Propagation,” Phys. Rev. D91, 053001 (2015), 
arXiv:1411.1174 [hep-ph] 

[34] Sergio Palomares-Ruiz, Aaron C. Vincent, and Olga 
Mena, “Spectral analysis of the high-energy Ice- 
Cube neutrinos,” Phys. Rev. D91, 103008 (2015), 
arXiv: 1502.02649 [astro-ph.HE] 

[35] M.G. Aartsen et al. (IceCube), “Elavor Ratio of As¬ 
trophysical Neutrinos above 35 TeV in IceCube,” 
Phys. Rev. Lett. 114, 171102 (2015), arXiv: 1502.03376 
[astro-ph.HE] 

[36] Andrea Palladino and Erancesco Vissani, “New 
parametrization of cosmic neutrino oscillations,” 
(2015), arXiv:1504.05238 [hep-ph] 

[37] M. G. Aartsen et aL (IceCube), “A combined maximum- 
likelihood analysis of the high-energy astrophysical neu¬ 
trino flux measured with IceCube,” Astrophys. J. 809, 
98 (2015), arXiv:1507.03991 [astro-ph.HE] 

[38] Walter Winter, “Describing the Observed Cosmic 
Neutrinos by Interactions of Nuclei with Matter,” 
Phys. Rev. D90, 103003 (2014), arXiv:1407.7536 [astro- 
ph.HE] 

[39] A. Palladino, G. Pagliaroli, E.L. Villante, and E. Vis¬ 
sani, “What is the Elavor of the Cosmic Neutrinos Seen 
by IceCube?” Phys. Rev. Lett. 114, 171101 (2015), 
arXiv: 1502.02923 [astro-ph.HE] 

[40] M.G. Aartsen et al. (IceCube), “IceCube-Gen2: A Vi¬ 
sion for the Euture of Neutrino Astronomy in Antarc¬ 
tica,” (2014), arXiv:1412.5106 [astro-ph.HE] 

[41] Raj Gandhi, Chris Quigg, Mary Hall Reno, and Ina 
Sarcevic, “Ultrahigh-energy neutrino interactions,” As- 
tropart. Phys. 5, 81-110 (1996), hep-ph/9512364 

[42] Raj Gandhi, Chris Quigg, Mary Hall Reno, and Ina 
Sarcevic, “Neutrino interactions at ultrahigh energies,” 
Phys. Rev. D58, 093009 (1998), hep-ph/9807264 

[43] Luis A. Anchordoqui, Haim Goldberg, Erancis Halzen, 
and Thomas J. Weiler, “Neutrinos as a diagnostic of 
high energy astrophysical processes,” Phys. Lett. B621, 
18-21 (2005), arXiv:hep-ph/0410003 [hep-ph] 

[44] Atri Bhattacharya, Raj Gandhi, Werner Rodejohann, 
and Atsushi Watanabe, “The Glashow resonance at Ice- 
Cube: signatures, event rates and pp vs. pj interac¬ 
tions,” JCAP 1110, 017 (2011), arXiv:1108.3163 [astro- 
ph.HE] 

[45] V. Barger, Lingjun Eu, J.G. Learned, D. Marfatia, 
S. Pakvasa, et al., “Glashow resonance as a window 
into cosmic neutrino sources,” Phys. Rev. D90, 121301 
(2014), arXiv: 1407.3255 [astro-ph.HE] 

[46] John G. Learned and Sandip Pakvasa, “Detecting tan- 


neutrino oscillations at PeV energies,” Astropart. Phys. 
3, 267-274 (1995), arXiv:hep-ph/9405296 

[47] John E. Beacom, Nicole E. Bell, Dan Hooper, Sandip 
Pakvasa, and Thomas J. Weiler, “Measuring flavor ra¬ 
tios of high-energy astrophysical neutrinos,” Phys. Rev. 
D68, 093005 (2003), erratum-ibid.D72, 019901 (2005), 
hep-ph/0307025 

[48] Edgar Bugaev, Teresa Montaruli, Yuri Shlepin, and 
Igor A. Sokalski, “Propagation of tan neutrinos and tan 
leptons through the earth and their detection in under¬ 
water / ice neutrino telescopes,” Astropart. Phys. 21, 
491-509 (2004), arXiv:hep-ph/0312295 [hep-ph] 

[49] Tamar Kashti and Eli Waxman, “Elavoring astrophys¬ 
ical neutrinos: Elavor ratios depend on energy,” Phys. 
Rev. Lett. 95, 181101 (2005), arXiv:astro-ph/0507599 

[50] M. Kachelriess, S. Ostapchenko, and R. Tomas, “High 
energy neutrino yields from astrophysical sources. 2. 
Magnetized sources,” Phys. Rev. D77, 023007 (2008), 
arXiv:0708.3047 [astro-ph] 

[51] S. Hummer, M. Maltoni, W. Winter, and C. Yaguna, 
“Energy dependent neutrino flavor ratios from cosmic 
accelerators on the Hillas plot,” Astropart. Phys. 34, 
205-224 (2010), arXiv: 1007.0006 [astro-ph.HE] 

[52] M. Kachelriess and R. Tomas, “High energy neutrino 
yields from astrophysical sources. 1. Weakly magnetized 
sources,” Phys. Rev. D74, 063009 (2006), arXiv:astro- 
ph/0606406 

[53] Rikard Enberg, Mary Hall Reno, and Ina Sarce¬ 
vic, “High energy neutrinos from charm in astro- 
physical sources,” Phys. Rev. D79, 053006 (2009), 
arXiv:0808.2807 [astro-ph] 

[54] S. Hummer, M. Riiger, E. Spanier, and W. Win¬ 
ter, “Simplified models for photohadronic interactions 
in cosmic accelerators,” Astrophys. J. 721, 630-652 
(2010), arXiv: 1002.1310 [astro-ph.HE] 

[55] C. Lunardini and A. Yu. Smirnov, “High-energy neu¬ 
trino conversion and the lepton asymmetry in the uni¬ 
verse,” Phys. Rev. D64, 073006 (2001), arXiv:hep- 
ph/0012056 [hep-ph] 

[56] Soebur Razzaque and A. Yu. Smirnov, “Elavor conver¬ 
sion of cosmic neutrinos from hidden jets,” JHEP 1003, 
031 (2010), arXiv:0912.4028 [hep-ph] 

[57] Sarira Sahu and Bing Zhang, “Effect of Resonant Neu¬ 
trino Oscillation on TeV Neutrino Elavor Ratio from 
Choked GRBs,” Res. Astron. Astrophys. 10, 943-949 
(2010), arXiv: 1007.4582 [hep-ph] 

[58] Yasaman Earzan and Alexei Yu. Smirnov, “Coherence 
and oscillations of cosmic neutrinos,” Nucl. Phys. B805, 
356-376 (2008), arXiv:0803.0495 [hep-ph] 

[59] Evgeny Akhmedov, Daniel Hernandez, and Alexei 

Smirnov, “Neutrino production coherence and os¬ 
cillation experiments,” JHEP 1204, 052 (2012), 

arXiv:1201.4128 [hep-ph] 

[60] B.J.P. Jones, “Dynamical pion collapse and the coher¬ 
ence of conventional neutrino beams,” Phys. Rev. D91, 
053002 (2015), arXiv: 1412.2264 [hep-ph] 

[61] K.A. Olive et al. (Particle Data Group), “Review of Par¬ 
ticle Physics,” Chin.Phys. C38, 090001 (2014) 

[62] M.C. Gonzalez-Garcia, Michele Maltoni, and Thomas 
Schwetz, “Updated fit to three neutrino mixing: sta¬ 
tus of leptonic CP violation,” JHEP 1411, 052 (2014), 
arXiv: 1409.5439 [hep-ph] 

[63] R.H. Dalitz, “On the analysis of tan-meson data and 
the nature of the tan-meson,” Phil. Mag. 44, 1068-1080 







7 


(1953) 

[64] Joachim Kopp, Pedro A. N. Machado, Michele Mal- 
toni, and Thomas Schwetz, “Sterile Neutrino Oscilla¬ 
tions: The Global Picture,” JHEP 1305, 050 (2013), 
arXiv:1303.3011 [hep-ph] 

[65] C. Giunti, M. Laveder, Y.F. Li, and H.W. Long, “Prag¬ 
matic View of Short-Baseline Neutrino Oscillations,” 
Phys. Rev. D88, 073008 (2013), arXiv:1308.5288 [hep- 
ph] 

[66] Sheldon L. Glashow, “Resonant Scattering of Antineu¬ 
trinos,” Phys. Rev. 118, 316-317 (1960) 

[67] V.S. Berezinsky and A.Z. Gazizov, “Gosmic Neutri¬ 
nos and Possibility to Search for W Bosons Having 
30-GeV-100-GeV Masses in Underwater Experiments,” 
JETP Lett. 25, 254-256 (1977). 

[68] V.S. Berezinsky and A.Z. Gazizov, “Neutrino - electron 
scattering at energies above the W boson production 
threshold,” Sov. J. Nucl. Phys. 33, 120-125 (1981). 

[69] H. Athar, G. Parente, and E. Zas, “Prospects for 
observations of high-energy cosmic tan neutrinos,” 
Phys. Rev. D62, 093010 (2000), arXiv:hep-ph/0006123 
[hep-ph] 

[70] Tyce DeYoung, S. Razzaque, and D.F. Gowen, 
“Astrophysical tan neutrino detection in kilometer- 
scale Gherenkov detectors via muonic tan decay,” 
Astropart. Phys. 27, 238-243 (2007), arXiviastro- 
ph/0608486 [astro-ph] 

[71] Zhi-zhong Xing and Shun Zhou, “The Glashow reso¬ 
nance as a discriminator of UHE cosmic neutrinos orig¬ 
inating from p-gamma and p-p collisions,” Phys. Rev. 
D84, 033006 (2011), arXiv:1105.4114 [hep-ph] 

[72] Atri Bhattacharya, Raj Gandhi, Werner Rodejohann, 
and Atsushi Watanabe, “On the interpretation of Ice- 
Gube cascade events in terms of the Glashow reso¬ 
nance,” (2012), arXiv:1209.2422 [hep-ph] 

[73] R. Abbasi et al. (IceGube), “A Search for UHE Tan Neu¬ 
trinos with IceGube,” Phys. Rev. D86, 022005 (2012), 
arXiv: 1202.4564 [astro-ph.HE] 

[74] V. Barger, J. Learned, and S. Pakvasa, “IceGube PeV 
cascade events initiated by electron-antineutrinos at 
Glashow resonance,” Phys. Rev. D87, 037302 (2013), 
arXiv: 1207.4571 [astro-ph.HE] 

[75] Mattias Blennow and Davide Meloni, “Non-standard in¬ 
teraction effects on astrophysical neutrino fluxes,” Phys. 
Rev. D80, 065009 (2009), arXiv:0901.2110 [hep-ph] 

[76] M. Bustamante, A.M. Gago, and Joel Jones Perez, 
“SUSY Renormalization Group Effects in Ultra 
High Energy Neutrinos,” JHEP 1105, 133 (2011), 
arXiv:1012.2728 [hep-ph] 

[77] Pankaj Jain, Supriya Kar, Douglas W. McKay, Sukanta 
Panda, and John P. Ralston, “Angular dependence 
of neutrino flux in KM^ detectors in low scale gravity 
models,” Phys. Rev. D66, 065018 (2002), arXivihep- 
ph/0205052 [hep-ph] 

[78] Shahid Hussain and Douglas W. McKay, “Energy and 
angular distribution of upward UHE neutrinos and sig¬ 
nals of low scale gravity: Role of tan decay,” Phys. Rev. 
D69, 085004 (2004), arXiv:hep-ph/0310091 [hep-ph] 

[79] E. Borriello, A. Guoco, G. Mangano, G. Miele, S. Pas¬ 
tor, et al., “Disentangling neutrino-nucleon cross sec¬ 
tion and high energy neutrino flux with a km^ neu¬ 
trino telescope,” Phys. Rev. D77, 045019 (2008), 
arXiv:0711.0152 [astro-ph] 

[80] D. Marfatia, D. W. McKay, and T. J. Weiler, “New 


physics with ultra-high-energy neutrinos,” Phys. Lett. 
B748, 113-116 (2015), arXiv: 1502.06337 [hep-ph] 

[81] John F. Beacom, Nicole F. Bell, Dan Hooper, Sandip 
Pakvasa, and Thomas J. Weiler, “Decay of high-energy 
astrophysical neutrinos,” Phys. Rev. Lett. 90, 181301 
(2003), arXiv:hep-ph/0211305 [hep-ph] 

[82] Michele Maltoni and Walter Winter, “Testing neutrino 
oscillations plus decay with neutrino telescopes,” JHEP 
07, 064 (2008), arXiv:0803.2050 [hep-ph] 

[83] Philipp Baerwald, Mauricio Bustamante, and Walter 
Winter, “Neutrino Decays over Gosmological Distances 
and the Implications for Neutrino Telescopes,” JGAP 
1210 , 020 (2012), arXiv: 1208.4600 [astro-ph.GO] 

[84] G. Pagliaroli, A. Palladino, F. Vissani, and F.L. Vil- 
lante, “Testing neutrino decay scenarios with IceGube 
data,” (2015), arXiv:1506.02624 [hep-ph] 

[85] Edward W. Kolb and Michael S. Turner, “Supernova 
SN 1987a and the Secret Interactions of Neutrinos,” 
Phys. Rev. D36, 2895 (1987) 

[86] Kunihto loka and Kohta Murase, “IceGube PeV-EeV 
neutrinos and secret interactions of neutrinos,” PTEP 
2014, 061E01 (2014), arXiv: 1404.2279 [astro-ph.HE] 

[87] Kenny G. Y. Ng and John F. Beacom, “Gosmic neutrino 
cascades from secret neutrino interactions,” Phys. Rev. 
D90, 065035 (2014), arXiv: 1404.2288 [astro-ph.HE] 

[88] Khr Blum, Anson Hook, and Kohta Murase, “High 
energy neutrino telescopes as a probe of the neutrino 
mass mechanism,” (2014), arXiv: 1408.3799 [hep-ph] 

[89] John F. Gherry, Alexander Friedland, and Ian M. Shoe¬ 
maker, “Neutrino Portal Dark Matter: From Dwarf 
Galaxies to IceGube,” (2014), arXiv: 1411.1071 [hep- 
ph] 

[90] Ayuki Kamada and Hai-Bo Yu, “Goherent Propagation 
of PeV Neutrinos and the Dip in the Neutrino Spectrum 
at IceGube,” (2015), arXiv:1504.00711 [hep-ph] 

[91] Anthony DiFranzo and Dan Hooper, “Searching for 
MeV-Scale Gauge Bosons with IceGube,” (2015), 
arXiv: 1507.03015 [hep-ph] 

[92] John F. Beacom et al, “Pseudo-Dirac neutrinos, a chal¬ 
lenge for neutrino telescopes,” Phys. Rev. Lett. 92, 
011101 (2004), arXiv:hep-ph/0307151 

[93] Arman Esmaili and Yasaman Farzan, “Implications 
of the Pseudo-Dirac Scenario for Ultra High En¬ 
ergy Neutrinos from GRBs,” JGAP 1212 , 014 (2012), 
arXiv:1208.6012 [hep-ph] 

[94] Anjan S. Joshipura, Subhendra Mohanty, and Sandip 
Pakvasa, “Pseudo-Dirac neutrinos via a mirror world 
and depletion of ultrahigh energy neutrinos,” Phys. Rev. 
D89, 033003 (2014), arXiv:1307.5712 [hep-ph] 

[95] John R. Ellis, J.S. Hagelin, Dimitri V. Nanopoulos, and 
M. Srednicki, “Search for Violations of Quantum Me¬ 
chanics,” Nucl. Phys. B241, 381 (1984) 

[96] Tom Banks, Leonard Susskind, and Michael E. Peskin, 
“Difficulties for the Evolution of Pure States Into Mixed 
States,” Nucl. Phys. B244, 125 (1984) 

[97] F. Benatti and R. Floreanini, “Open system approach to 
neutrino oscillations,” JHEP 02, 032 (2000), arXiv:hep- 
ph/0002221 

[98] A.M. Gago, E.M. Santos, W.J.G. Teves, and 
R. Zukanovich Funchal, “A Study on quantum decoher¬ 
ence phenomena with three generations of neutrinos,” 
(2002), arXiv:hep-ph/0208166 [hep-ph] 

[99] Dean Morgan, Elizabeth Winstanley, Jurgen Brunner, 
and Lee F. Thompson, “Probing quantum decoher- 








ence in atmospheric neutrino oscillations with a neu¬ 
trino telescope,” Astropart. Phys. 25, 311-327 (2006), 
arXiv:astro-ph/0412618 

[100] Dan Hooper, Dean Morgan, and Elizabeth Winstanley, 
“Lorentz and CPT invariance violation in high-energy 
neutrinos,” Phys. Rev. D72, 065009 (2005), arXivihep- 
ph/0506091 

[101] Luis A. Anchordoqui, Haim Goldberg, M.C. Gonzalez- 
Garcia, Francis Halzen, Dan Hooper, et aL, “Probing 
Planck scale physics with IceGube,” Phys. Rev. D72, 
065019 (2005), arXiv:hep-ph/0506168 [hep-ph] 

[102] Mauricio Bustamante, John F. Beacom, and Kohta 
Murase, In preparation. 

[103] Don Golladay and V. Alan Kostelecky, “Lorentz violat¬ 
ing extension of the standard model,” Phys. Rev. D58, 
116002 (1998), arXiv:hep-ph/9809521 [hep-ph] 

[104] V. Alan Kostelecky and Matthew Mewes, “Lorentz and 
GPT violation in the neutrino sector,” Phys. Rev. D70, 
031902 (2004), arXiv:hep-ph/0308300 [hep-ph] 

[105] John Ellis and Nick E. Mavromatos, “Probes of 
Lorentz Violation,” Astropart. Phys. 43, 50-55 (2013), 
arXiv:1111.1178 [astro-ph.HE] 

[106] M. Bustamante, A.M. Gago, and G. Pena-Garay, 
“Energy-independent new physics in the flavour ratios 
of high-energy astrophysical neutrinos,” JHEP 1004, 
066 (2010), arXiv: 1001.4878 [hep-ph] 

[107] M. Gasperini, “Experimental Gonstraints on a Minimal 
and Nonminimal Violation of the Equivalence Principle 
in the Oscillations of Massive Neutrinos,” Phys. Rev. 
D39, 3606-3611 (1989) 

[108] M.N. Butler, S. Nozawa, R.A. Malaney, and A.I. 
Boothroyd, “Gravitationally induced neutrino oscilla¬ 
tions,” Phys. Rev. D47, 2615-2618 (1993) 

[109] S.L. Glashow, A. Halprin, P.I. Krastev, Ghung Ngoc Le¬ 
ung, and James T. Pantaleone, “Gomments on neutrino 
tests of special relativity,” Phys. Rev. D56, 2433-2434 
(1997), arXiv:hep-ph/9703454 [hep-ph] 

[110] V. De Sabbata and M. Gasperini, “Neutrino Oscillations 
in the Presence of Torsion,” Nuovo Gim. A65, 479-500 
(1981) 

[111] Garlos A. Argelles, Teppei Katori, and Jordi Sal- 
vado, “New Physics in Astrophysical Neutrino Flavor,” 
(2015), arXiv:1506.02043 [hep-ph] 

[112] Elke Aeikens, Heinrich Pas, Sandip Pakvasa, and 
Philipp Sicking, “Flavor ratios of extragalactical neu¬ 
trinos and neutrino shortcuts in extra dimensions,” 
(2014), arXiv:1410.0408 [hep-ph] 

[113] Kaustubh Agashe, Yanou Gui, Lina Necib, and Jesse 
Thaler, “(In)direct Detection of Boosted Dark Matter,” 
JGAP 1410, 062 (2014), arXiv: 1405.7370 [hep-ph] 

[114] Joachim Kopp, Jia Liu, and Xiao-Ping Wang, “Boosted 
Dark Matter in IceGube and at the Galactic Genter,” 
JHEP 1504, 105 (2015), arXiv: 1503.02669 [hep-ph] 

[115] Poonam Mehta and Walter Winter, “Interplay of en¬ 
ergy dependent astrophysical neutrino flavor ratios 
and new physics effects,” JGAP 1103, 041 (2011), 
arXiv:1101.2673 [hep-ph] 

[116] A. Kappes (KM3NeT), “KM3NeT: A Next Gen¬ 
eration Neutrino Telescope in the Mediterranean 
Sea,” Proceedings of the Sixth International Work¬ 
shop on New Worlds in Astroparticle Physics^ (2007), 
arXiv:0711.0563 [astro-ph] 

[117] Svenja Hummer, Philipp Baerwald, and Walter Win¬ 
ter, “Neutrino Emission from Gamma-Ray Burst Fire¬ 


balls, Revised,” Phys. Rev. Lett. 108, 231101 (2012), 
arXiv:1112.1076 [astro-ph.HE] 

[118] P. D. Serpico and M. Kachelriess, “Measuring the 13- 
mixing angle and the GP phase with neutrino tele¬ 
scopes,” Phys. Rev. Lett. 94, 211102 (2005), hep- 
ph/0502088 

[119] Walter Winter, “How astrophysical neutrino sources 
could be used for early measurements of neutrino mass 
hierarchy and leptonic GP phase,” Phys. Rev. D74, 
033015 (2006), hep-ph/0604191 

[120] Kfir Blum, Yosef Nir, and Eli Waxman, “Probing 
GP violation in neutrino oscillations with neutrino tele¬ 
scopes,” (2007), arXiv:0706.2070 [hep-ph] 

[121] Davide Melon! and Tommy Ohlsson, “Leptonic GP vi¬ 
olation and mixing patterns at neutrino telescopes,” 
Phys. Rev. D86, 067701 (2012), arXiv:1206.6886 [hep- 
ph] 

[122] Animesh Ghatterjee, Moon Moon Devi, Monojit Ghosh, 
Reetanjali Moharana, and Sushant K. Rant, “Probing 
GP violation with the first three years of ultrahigh en¬ 
ergy neutrinos from IceGube,” Phys. Rev. D90, 073003 
(2014), arXiv:1312.6593 [hep-ph] 

[123] Lingjun Fu and Ghiu Man Ho, “Flavor Ratios and 
Mass Hierarchy at Neutrino Telescopes,” (2014), 
arXiv: 1407.1090 [hep-ph] 

[124] G. Adams et al (LBNE), “The Long-Baseline Neutrino 
Experiment: Exploring Fundamental Symmetries of the 
Universe,” (2013), arXiv:1307.7335 [hep-ex] 

[125] John F. Beacom, Nicole F. Bell, Dan Hooper, Sandip 
Pakvasa, and Thomas J. Weiler, “Sensitivity to O 13 
and S in the decaying astrophysical neutrino scenario,” 
Phys. Rev. D69, 017303 (2004), arXiv:hep-ph/0309267 
[hep-ph] 









9 


Supplemental Material 


Appendix A: Plots for inverted mass hierarchy 


We repeat the figures shown in the main text of the 
paper, but now for the inverted hierarchy (IH) instead 
of the normal hierarchy. The differences are modest and 
are vanishing for the 3a regions. 



FIG. 5. Same as Fig.0 but for an inverted mass hierarchy. FIG. 7. Same as Fig. but for an inverted mass hierarchy. 
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FIG. 6. Same as Fig. but for an inverted mass hierarchy. FIG. 8. Same as Fig.[^ but for an inverted mass hierarchy. 
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Appendix B: Energy-dependent effects on the flavor 
composition 

The energy dependence of the flavor composition at 
the source has been considered in, e.g., Refs. j49ti5^ . 
At lower energies, a | : 0 )g composition is expected, 
coming from the full pion decay chain. At higher en¬ 
ergies, where the synchrotron losses of the muons cre¬ 
ated by the pion decays become larger and where addi¬ 
tional neutrino production channels become accessible, 
the source flavor composition changes. Non-trivial flavor 
compositions are obtained if several processes compete 
or the cooled muons pile up at lower energies m- 

If neutrinos in a broad energy interval are considered, 
then the inferred flavor composition at Earth will be a 
superposition of those at different energies. It will be 
challenging for IceCube to pinpoint the exact flavor ratios 
unless a large volume of data is available for fine energy 
binning to be feasible. 

Let us assume that such a binning is in fact feasible. 
In that case, by assuming that the energy dependence of 
the source flavor composition is the same for all of the 
sources contributing to the neutrino diffuse flux, we can 
predict the ratios at Earth at different energies. 

Eigure shows the variation of the flavor ratios at 
Earth from a source with emission parameters given by 
a “classical” pion beam source evolving into a muon- 
damped source at higher energies (test point TP 13 from 
Ref. [51], a photohadronic model where the target pho¬ 
tons come from synchrotron emission of co-accelerated 
electrons). As expected, the trajectory in flavor space 
lies within the Standard Model region shown in Eig. [^ 

Eigure shows the results for a different parameter 
set (test point TPS from Ref. [51] , corresponding to AGN 
cores, a mixed source at low energies, where muon pile- 
ups add to the pion decay chain, evolving into muon 
damped source at high energies), this time taking into 
account in addition invisible decays of U 2 and z/ 3 , and a 
stable z/i, z.e., in the normal hierarchy. The fact that 


decays are invisible means that the decay products of 1^2 
and z /3 do not contribute to the flux. 


In this example, the energy dependence of the fla¬ 
vor composition at the source competes with the energy 
dependence of the new physics effect; see discussion in 
Ref. mg. We have fixed the decay damping parameter 
to aL = 10^ GeV, following the notation in Eq. (18) of 
Ref. |115j . As expected, the trajectory in flavor space 
now leaves the Standard Model region at the lowest en¬ 
ergies, where the suppression from decay is stronger, and 
reaches the region corresponding to the flavor content of 
ui (see Eig. Q. The full trajectory is still completely 
contained within the new physics region from Eig. 


Eigure 11 shows complementary results for test point 
TPS in the inverted hierarchy case, when z /3 is stable 
while and U 2 decay. In this case, the trajectory leaves 
the Standard Model region and reaches the region corre¬ 
sponding to the flavor content of Note that this case 
may be disfavored by the observation of neutrinos from 
SN 1987A. 


As long as the flavor ratios at Earth are confined within 
the gray region of Eigs. [qHTT] no distinction between 
the effects of energy-dependent new physics and those 
of energy-dependent flavor ratios at the source is pos¬ 
sible. Given sufficient energy resolution and statistics, 
however, it might become possible to disentangle them: 
if, for some energy range, the flavor ratios leave the gray 
region and enter the region that is reachable exclusively 
with new physics, shown in Eigs.|^and[^ then this would 
represent evidence of energy-dependent new physics. As 
neutrino mixing parameters become better measured, the 
gray region will shrink, improving the ability to deter¬ 
mine the cause of any observed energy variation. 

If the experimental determination of the flavor ratios 
finds no energy dependence, then we could interpret this 
as a hint that neutrino production occurs in a narrow 
energy window, inside of which the flavor composition at 
the sources is approximately constant. 
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FIG. 9. Left: Variation with energy of the flavor ratios at Earth, for neutrinos coming from a source with emission parameters 
given by test point TP 13 from Ref. m- The neutrinos were produced in py interactions, calculated using NeuCosmA [5^1117] . 
Right: Ternary plots showing the trajectory of the flavor ratios of the same source, marking the values at different energies. 



FIG. 10. Left: Variation with energy of the flavor ratios at Earth, for neutrinos coming from a source with emission parameters 
given by test point TPS from Ref. m, including decay of 1^2 and 1 ^ 3 . The neutrinos were produced in py interactions, calculated 
using NeuGosmA [5^ 1117] . Right: Ternary plots showing the trajectory of the flavor ratios of the same source, marking the 
values at different energies. 
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FIG. 11. Left: Variation with energy of the flavor ratios at Earth, for neutrinos coming from a source with emission parameters 
given by test point TPS from Ref. m, including decay of ui and 1^2 ■ The neutrinos were produced in py interactions, calculated 
using NeuCosmA [54l 1117] . Right: Ternary plots showing the trajectory of the flavor ratios of the same source, marking the 
values at different energies. 
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FIG. 12. Flavor ratios at Earth with unconstrained source flavor ratios and four particular choices of source flavor ratios. The 
mixing angles are kept fixed at their current best-fit values; the CP-violation phase (5cp is the only mixing parameter that has 
been varied to create the colored regions. The regions corresponding to the flavor content of v\ and at the la and 3a regions 
are also shown; note that the flavor content of has no dependence on (5cp. 


Appendix C: Measurement of fcp? 

There have been proposals to determine the leptonic 
CP phase, (5cp, using high-energy astrophysical neutri¬ 
nos; see, e.^.. Refs. |11814123] . These studies considered 
specific flavor ratios at the source, typically, (| : | : 0)g, 
(0 : 1 : 0)g, or (1 : 0 : 0)g, and studied the effect of ^cp on 
the ratio at Earth of the flux to the 

flux. For a value of Ois ~ 10°, this ratio was found to 
vary as a function of ^cp in np to ^ 20%, with the muon- 
damped and neutron decay scenarios being the most sen¬ 
sitive. At higher energies, if tan-neutrinos are identified, 
the ratio of + Pg lo can be used to improve 

the measurement of ^cp- While optimistic studies inni 
find that astrophysical sources by themselves might be 
able to claim a measurement of (5cp, most studies (e.^., 
|118[I119] 1 focus on a combined analysis of oscillation and 
astrophysical neutrino data in order to ascertain discov¬ 
ery. 

We therefore show in Fig. the allowed regions of 
flavor composition if all of the mixing parameters except 
^cp are fixed to their current best-fit values. This as¬ 
sumption corresponds to future improved measurements 
of the mixing angles coming from other experiments. 


whereas Sqp is to be determined from neutrino telescope 
data. We include the region corresponding to uncon¬ 
strained flavor ratios at the source, not only to particular 
ratios. For comparison, we include the outlines of the re¬ 
gions generated when the mixing angles are also allowed 
to vary within their la and 3a regions. 

The comparison of Figs. [^&[^and Fig, shows that, 
in the absence of knowledge of the flavor composition 
at the source, 6 cp cannot be extracted. However, if 
the source flavor composition were known to be muon- 
damped ((0:1: 0)s) or neutron source ((1:0: 0)s), 
then, given the statistics of IceCube-Gen2, one would 
expect a marginal statistically significant contribution to 
the mass hierarchy or CP-violation sensitivities of T2K 
and NOpA m- A new long-baseline experiment, such 
as DUNE |124j . would be more precise. 

Notably, if neutrinos fully decay in the NH, the mea¬ 
surement of ^cp iiiay even be more attractive [US]. Fig¬ 
ure!^ shows that, in this case, IceCube-Gen2 could actu¬ 
ally measure the phase, since the flavor ratios would be 
determined solely by the flavor content of the one stable 
mass eigenstate, ui (lower right yellow region). This de¬ 
cay scenario is, however, disfavored by current IceCube 
data at > 2a; see Fig. 


























